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Abstract 

In this work we compare two different approaches to calculation of the three-body 
resonances on the basis of Faddeev differential equations. The first one is the com- 
plex scaling approach. The second method is based on an immediate calculation of 
resonances as zeros of the three-body scattering matrix continued to the physical 
sheet. 
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1 Introduction 



The complex scaling method [1,2] invented in early 70-s remains one of the 
most effective approaches to calculation of resonances in few-body systems. 
This method is applicable to an iV-body problem in the case where inter- 
particle interaction potentials are analytic functions of coordinates. The com- 
plex scaling gives a possibility to rotate the continuous spectrum of an iV-body 
Hamiltonian in such a way that certain sectors of unphysical sheets neighbor- 
ing the physical one turn into a part of the physical sheet for the resulting 
non- self adjoint operator. Resonances appear to be complex eigenvalues of this 
operator [1,2] while the binding energies stay fixed during the scaling trans- 
formations. Therefore, when searching for the resonances within the complex 
scaling approach one may apply the methods which are usually employed to 
locate the binding energies. Some reviews of the literature on the complex 
scaling and its many applications can be found, in particular, in [4-6]. Here 
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we only mention that there is a rigorous mathematical proof [3] that for a 
rather wide class of interaction potentials the resonances given by the com- 
plex scaling method coincide with the "true scattering resonances" , i. e. the 
poles of the analytically continued scattering matrix in the unphysical sheets. 

Along with the complex scaling, various different methods are also used for 
calculations of the resonances. Among the methods developed to calculate 
directly the scattering-matrix resonances we, first, mention the approach based 
on the momentum space Faddeev integral equations [7,8] (see, e. g., Ref. [9] and 
references cited therein). In this approach one numerically solves the equations 
continued into an unphysical sheet and, thus, the three-body resonances arise 
as the poles of the continued T-matrix. Another approach to calculation of the 
scattering-matrix resonances is based on the explicit representations [10,11] for 
the analytically continued T- and S-matrices in terms of the physical sheet. 
From these representations one infers that the three-body resonances can be 
found as zeros of certain truncations of the scattering matrix only taken in 
the physical sheet. Such an approach can be employed even in the coordinate 
space [11,12]. 

To the best of our knowledge there are no published works applying the com- 
plex scaling to the Faddeev equations. Therefore, we consider the present 
investigation as a first attempt undertaken in this direction. However, the pur- 
pose of our work is rather two-fold. On the one hand, we make the complex 
scaling of the Faddeev differential equations. On the other hand we compare 
the complex scaling method with the scattering-matrix approach suggested 
in [11,12]. We do this making use of both the approaches to examine reso- 
nances in a model system of three bosons having the nucleon masses and in 
the three-nucleon (nnp) system itself. 



2 Formalism 

First, we recall that, after the scaling transformation, the three-body Schrodinger 
operator reads as follows [1-3] 



where $ = iO is the scaling parameter with 6 G M. By Ax we understand 
the six-dimensional Laplacian in X = (x a ,y a ) where x a ,y a are the standard 
Jacobi variables, a — 1, 2, 3. Notation V a is used for the two-body potentials 
which are assumed to depend on x a = |x Q | but not on x a = x a /x a . 

The corresponding scaled Faddeev equations which we solve read 
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(1) 



a=l 



2 



[-e-™A x + v a (e»x a ) - z]& a \z;X) + V a (e*x a ) £ ^ W) (z;X) = f a (X), 

a = 1,2, 3. (2) 

Here / = (/i, / 2 , /a) is an arbitrary three-component vector with components 
/ a belonging to the three-body Hilbert space L 2 (M 6 ). 

The partial-wave version of the equations (2) for a system of three identical 
bosons at the zero total angular momentum L = reads 

e- 2ie H^^( Z] x,y)-z x, y) + V(e*a;)tt t (*; x, y) = f (l \x, y) (3) 
where x > 0, y > and if^ denotes the partial- wave kinetic energy operator, 



& = -&-W +,{, + l) ? + ^)- ' = 0.2.4,..., 



while stands for the partial-wave component of the total wave function, 

+1 

x, y) = $ t (z] x,y) / drj h w (x, y, rj) $ v (z; x', y') . (4) 

v _i 



Here, x' = y | x 2 + | y 2 — ^ xy?? and y' = y | x 2 + | y 2 + xy??. Explicit 
expression for the geometric function hw(x,y,r]) can be found, e.g., in [8]. 

The partial-wave equations (3) are supplied with the boundary conditions 

&i(z;x,y)\ x=o = and x, y)\ y=0 = 0. (5) 

For compactly supported inhomogeneous terms /^(x, y) the partial- wave Fad- 
deev component $i(z;x, y) also satisfies the asymptotic condition 



$i(z; x, y) = 5i ipd(e ie x) exp(iV^ - e d e l9 y) £l (z) + o (y 1/2 ) 

'A l ( Z ;y/x) + o(p- 1 / 2 ) 



(6) 



For simplicity it is assumed in this formula that the two-boson subsystem 
has only one bound state with the energy e^, and ipd( x ) represents its wave 
function. The values of ao and Ai(y/x) are the main asymptotical coefficients 
effectively describing the contributions to from the elastic (2 + 1-^2 + 1) 
and breakup (2 + 1 — > 1 + 1 + 1) channels, respectively. Hereafter, by ^/C 
( G C, we understand the main (arithmetic) branch of the function C 1 / 2 . 
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In the scaling method a resonance is looked for as the energy z which produces 
a pole to the quadratic form 

Q(6,z) = ([H F (6)-z}- 1 fj) 

where Hp (9) is the non-selfadjoint operator resulting from the complex-scaling 
transformation of the Faddeev operator. The latter operator is just the opera- 
tor constituted by the 1. h. s. parts of Eqs. (2). The resonance energies should 
not, of course, depend on the scaling parameter 9 and on the choice of the 
terms f^(x, y). 

In the scattering-matrix approach we solve the same partial-wave Faddeev 
equations (3) with the same boundary conditions (5) and (6) but for 9 = 
and 

+i 

fV\x,y) = -V{x) J drjhio(x : y : rj) tp d (x') sin(^z - e d y'). 

-i 

The resonances are looked for as zeroes of the truncated scattering-matrix (see 
[12] for details) sq(z) = l + 2iao(-2), where the (1 + 1 — > 1 + 1) elastic scattering 
amplitude &o(z) for complex energies z in the physical sheet is extracted from 
the asymptotics (6). 

For numerical solution of the boundary- value problem (3-6) we employ its 
finite-difference approximation in the hyperradius-hyperangle coordinates. A 
detail description of the finite-difference algorithm used can be found in Ref. 
[13]. 



3 Results 

In the table we present our results obtained for a complex-scaling resonance 
in the model three-body system which consists of identical bosons having the 
nucleon mass. To describe interaction between them we employ a Gauss-type 
potential of Ref. [12] 

V(x) = V exp[-fi x 2 ] + V b exp[-(x b (x - x b ) 2 ] 

with V = -55MeV, /i = 0.2 fm" 2 , x b = 5 fin, /i b = 0.01 fm" 2 and V b = 
1.5. The figures in the table correspond to the roots of the inverse function 
[Q(9, z)]~ l for L = and / = only taken into account. In the present 
calculation we have taken up to 400 knots in both hyperradius and hyper- 
angle variables while for the cut-off hyperradius we take 40 fm. One observes 
from the table that the position of the resonance depends very weakly on 
the scaling parameter 9 which confirms a good numerical quality of our re- 
sults. We compare the resonance values of the table to the resonance value 
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Fig. 1. Trajectory of the resonance z res in the model system of three bosons with the 
nucleon masses. Values of the barrier Vb in MeV are given near the points marked 
on the curve. 

^res = —5.952 — 0.403 i MeV obtained for the same three-boson system with 
exactly the same potentials but in the completely different scattering-matrix 
approach of Ref. [12]. We see that, indeed, both the complex scaling and the 
scattering matrix approaches give the same result. 



e 


z res (MeV) 


9 


z res (MeV) 


0.25 


-5.9525 - 0.4034 i 


0.50 


-5.9526 - 0.4032 i 


0.30 


-5.9526 - 0.4033 i 


0.60 


-5.9526 - 0.4033 i 


0.40 


-5.9526 - 0.4032 i 


0.70 


-5.9526 - 0.4034 i 



We also watched the trajectory of the above resonance when the barrier ampli- 
tude Vb varied (see. Fig. 1). While the complex scaling method was applicable 
it gave practically the same positions for the resonance. For the barrier ampli- 
tudes Vb smaller than 1.0 only the scattering-matrix approach allows to locate 
the resonance (which finally, for Vf, < 0.85, turns into a virtual level). 

As to the nnp system in the S- state where we employed the MTI-III [14] 
potential model, both the methods applied give no resonances on the two- 
body unphysical sheet (see [12]). Moreover, we have found no resonances in 
the part of the three-body sheet accessible via the complex scaling method. 
Thus, at least in the framework of the MT I — III model we can not confirm the 
experimental result of Ref. [16] in which the point — 1.5±0.3— i(0.3±0.15) MeV 
was interpreted as a resonance corresponding to an exited state of the triton 
3 H. 

The triton virtual state can be only calculated within the scattering-matrix 
method but not in the scaling approach. Our present improved scattering- 
matrix result for the triton virtual state is —2.690 MeV (i.e. the virtual level 
lies 0.47 MeV below the two-body threshold). This result has been obtained 
with the MT I-III potential on a grid having 1000 knots in both hyperradial 
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and hyperradial variables and with the value of cut-off hyperradius equal to 
f20fm. Notice that some values for the virtual-state energy obtained by dif- 
ferent authors can be found in [9] and all of these values are about 0.5 MeV 
below the two-body threshold. 
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